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Effects  of  percolation  on  the  energy  transfer  from  an  admoLecule  to  a 
two-component  granular  composite  substrate  are  investigated.  The  substrate 
material  is  a  random  mixture  of  a  metal  and  an  insulator  non- absorptive  to  the 
light.  The  luminescence  decay  rate  of  the  admolecule  near  the  composite 
surface  is  calculated  by  means  of  the  effective  medium  theory  as  a  function  of 
the  volume  fraction  of  the  metallic  granules.  It  is  found  that  the  decay  rate 
is  greatly  enhanced  around  the  dc  conduction  percolation  threshold  at  low 
transition  frequencies.  For  high  frequencies,  the  decay  rate  peaks  at  a 
higher  volume  fraction  of  the  metallic  particles,  indicating  the  existence  of 
a  "bulk  percolation"  threshold. 
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I.  Introduction 

It  is  well  known  that  in  a  two -component  particle  composite  material, 

the  conductivity  exhibits  singular  behavior  near  a  percolation  threshold, 

namely,  a  composition  at  which  one  of  the  constituents  forms  a  connected  path 

extending  throughout  the  material.^  As  an  example,  the  effective  conductivity 

of  a  composite  material  composed  of  metal  and  insulator  particles  varies  as 

near  the  metal  percolation  threshold  ,  where  q  is  the  volume 

fraction  of  the  metal  and  t  is  a  characteristic  exponent  which  remains  to  be 

determined.  Similar  behavior  has  also  been  found  for  the  nonlinear 

2 

susceptibility  in  the  static  limit. 

In  addition  to  static  properties,  optical  properties  of  the  two- 

component  composite  material,  both  linear  and  nonlinear,  have  been  extensively 
3-9 

studied.  In  general,  composite  materials  have  optical  properties  very 
different  from  those  of  their  constituents,  especially  in  the  vicinity  of  the 

percolation  threshold.  For  instance,  the  observation  of  a  large  enhancement 

6  8  7 

(10  -10  )  of  the  effective  optical  nonlinearity  has  been  reported  in 

composite  materials.  It  is  therefore  interesting  to  investigate  their  optical 

properties  near  the  percolation  threshold. 

On  the  other  hand,  optical  properties  of  an  atomic  system  located  near  a 

solid  surface  can  change  qualitatively  depending  upon  the  structure  of  the 

substrate.  Detailed  information  about  the  spectroscopic  properties  of  the 

adspecies  can  be  employed  to  probe  the  electronic  and  optical  properties  of 

the  substrate.  A  recent  study  of  the  decay  rate  and  resonance  fluorescence 

spectrum  of  a  molecule  adsorbed  near  the  surface  of  a  particle  composite^^ 

illustrates  clearly  the  strong  dependence  of  the  energy  transfer  on  the  volume 

fraction  of  the  impurity. 
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We  investigate  in  this  communication  the  effect  of  percolation  on  the 
energy  transfer  between  an  admolecule  and  the  substrate.  We  consider  a  two- 


component  granular  material,  that  is,  a  random  mixture  of  two  kinds  of  small 
particles;  one  is  an  insulator  not  absorptive  to  the  light,  and  the  other  is 
a  metal.  Optical  properties  of  this  material  are  described  by  an  effective 
dielectric  function  derived  by  means  of  the  effective  medium  theory 
(EMT).^’^^  When  one  component  of  the  composite  just  forms  an  infinite 
connected  path  through  the  material,  it  is  at  the  percolation  threshold. 
Thus,  the  conductivity  increases  suddenly  when  the  metal  granules  just  reach 
the  percolation  threshold.  Our  purpose  is  to  investigate  the  luminescence 
decay  rate  of  an  admolecule  near  the  surface  of  the  composite  which  is  in  the 
vicinity  of  its  percolation  threshold.  It  is  found  that  the  decay  rate  shows 
a  peak  centered  around  the  conduction  percolation  threshold  at  low 


frequencies,  or  when  u  «  where  is  the  plasmon  frequency  of  the  metal. 
For  adraolecules  of  high  transition  frequency,  or  w  -  the  peak  occurs  at  a 
higher  q,  indicating  the  possible  existence  of  another  percolation  threshold. 


Implications  of  this  and  other  results  will  be  discussed. 


II .  Theory 

We  consider  a  molecule  adsorbed  at  a  distance  d  from  the  surface  of  a 
two-component  particle  composite  material  as  the  substrate.  For  simplicity, 
we  assume  that  the  molecule  has  only  two  levels  |+  and  | ->  with  an  energy 
difference  -  E  - 

The  optical  spectrum  of  the  scattered  field  from  an  admolecule  can  in 

general  be  treated  by  means  of  a  set  of  surface-dressed  optical  Bloch 
14 

equations  Including  effects  of  the  surface-reflected  field.  For  the  present 
purpose  of  the  luminescence  decay-rate  study,  it  is  sufficient  to  consider  the 


k 


admolecule  as  a  damped  harmonic  oscillator.  Thus,  the  Hamiltonian  of  the 
system  is  given  by 

H  -  (w-i-y  )a^a  -  p*a^E  ,  (1) 

O  IV 

where  a^(a)  is  the  creation  (annihilation)  operator  for  the  molecular  dipole 
oscillator  whose  matrix  element  is  denoted  by  p,  and  denotes  the  reflected 
electric  field  in  the  direction  of  the  dipole  moment.  Following  Dekker's 
quantization  procedure  for  dissipative  systems,  the  equation  of  motion  for 
the  density  operator  of  the  molecular  system  can  be  written  as 

;  -  >[a^[a,H]p]  +  ^[p[H^a^l.a]  .  (2) 

With  the  help  of  Eqs.  (1)  and  (2),  we  find  that  the  mean  value  of  the  operator 
obeys  the  equation 

~<a>  -  Tr<pa>  -  - i(c<H-A«)<a>  -  j<a>  ,  (3) 

at 

with  the  decay  rate  given  by 

2y^/co 

In  the  above  derivation,  we  have  made  use  of  the  following  definitions.  Both 
the  frequency  shift  Aw  and  width  change  Ay  -  7*7q  are  caused  by  the  presence 
of  the  surface.  The  decay  rate  of  the  excited  molecule  inside  a  homogeneous 
meditun  is  ^  "  /c  •  and  the  dipole  moment  operator  is  -  pa.  In 
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addition,  we  have  assumed  that  the  molecule  is  in  a  medium  of  dielectric 
constant  £^,  and  Kq  -  w/c. 

As  is  well  known,  the  reflected  field  for  an  emitting  dipole 
with  moment  located  near  the  surface  of  a  solid  substrate  can  be  found  by 
solving  Maxwell  equations  with  corresponding  boundary  conditions.  Since  the 
substrate  we  are  considering  now  is  a  two-component  particle  composite 
material  with  symmetric  microgeometry,  E_  must  depend  on  the  optical 
properties  of  both  components  as  well  as  their  volume  fractions.  Therefore 
measurements  of  the  decay  rate  alone  can  reveal  much  of  the  characteristic 
properties  of  the  substrate  via  E  .  Here  we  just  consider  the  case  in  which 
the  admolecule  is  oriented  such  that  its  dipole  moment  is  normal  to  the 
substrate  surface.  The  decay  rate  can  then  be  worked  out  from  Eq.  (4)  by  the 
standard  procedure, namely, 


q  +  (1-q) 


(6) 
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where  e  (e  )  stands  for  Che  dielectric  function  of  the  insulator  (metal),  q  is 
di  D 

the  volume  fraction  of  the  conductor  in  the  composite,  and  is  the  critical 
volume  fraction  of  percolation  for  the  connection  of  metal  particles.^  In  the 
standard  EMT,  q^  -  1/3. 

Equation  (6)  can  be  rewritten  as 


<q  -»'2  ■  (■’-V'bl'a 


€  e,  -  0 
a  b 


The  quadratic  equation  (7)  has  in  general  two  solutions  that  are  complex 

conjugate  to  each  other.  We  are  interested  in  a  physical  solution  £2^*^)  which 

is  a  continuous  function  of  u>  and  has  a  positive  imaginary  part.  Since  the 

percolation  threshold  of  a  metal- insulator  composite  material  changes  with  the 

microgeometry  structure  of  particles,  q  may  be  regarded  as  a  parameter 

o 

characterizing  the  microstructure  of  the  two  components. 
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typical  value  which  can  always  be  satisfied  by  adjusting  the  molecular 
distance  from  the  surface  for  a  given  transition  frequency  w. 

We  first  solve  Eq.  (7)  for  the  effective  dielectric  function  of  the 
composite  substrate  with  given  q^  and  w.  Plotted  as  functions  of  the  volume 
fraction  q,  both  the  real  and  imaginary  parts  of  shown  in  Fig.  1  for  q^ 
-0.3  and  three  typical  frequencies;  (a)  w  -  0.1  in  the  low-frequency 
region:  The  absorption  part  peaks  around  q  -  q^ ,  where  the  dispersion  part 
changes  its  sign.  Hence  the  substrate  is  metal- like  for  q  >  q^  in  low- 
frequency  regime.  Our  numerical  study  shows  similar  results  for  different 
values  of  q^ ;  (b)  u  -  0.35  :  Both  the  peak  of  the  imaginary  part  and  the 
zero  of  the  real  part  move  toward  larger  q;  and  (c)  w  -  :  The  real  part 
remains  positive  in  the  whole  range,  and  the  imaginary  part  has  a  small  peak 
around  q  »  0.7.  We  believe  that  this  is  due  to  the  bulk  effect  on  the 
percolation  threshold,  which  we  shall  discuss  later. 

Since  the  surface  plasmon  contribution  leads  to  a  minimum  in  the 
denominator  in  Eq.  (7),^^  this  resonance  condition  determines  the  dispersion 
relation  of  the  surface  plasmon.  If  the  substrate  is  nonabsorptive ,  the 
surface  plasmons  have  infinite  lifetime  and  £^^2  ^2^1  ”  of  the 
nonzero  in  the  present  case,  the  lifetime  remains  finite.  In  the  low- 
frequency  region,  say,  <  0.2,  the  absorption  part  of  the  optical 
susceptibility  of  the  composite  material  is  usually  very  large,  as  seen  in 
Fig.  1(a),  and  hence  the  damping  rate  is  too  large  for  surface  plasmons  to 
propagate.  In  the  high-frequency  region,  w/Wp  >  0.9,  it  is  seen  from  Fig. 
1(c)  that  the  real  part  of  «2  is  always  positive  even  though  the  imaginary 
part  is  much  smaller.  There  can  not  be  surface  plasmon  waves  here  either. 


because  the  real  part  of  £2  ®'ist  be  negative  (at  least)  for  surface  modes  to 
exist. 
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The  decay  rate  is  computed  from  (5)  with  £2(w)  given  by  (7).  Since  we 
are  interested  in  percolation  effects  on  the  energy  transfer  from  the 
admolecule  to  the  surface,  7  is  calculated  as  a  function  of  the  volume 
fraction  q  of  the  metal  particles  and  is  expressed  in  the  unit  of  7^.  It  is 
well  known  that  the  conduction  percolation  threshold  depends  sensitively 
upon  the  microstructure  of  the  two-component  system.  For  a  metal- insulator 
system  of  symmetric  microgeometry,  q^  is  not  expected  to  exceed  0.5  normally. 
Some  of  our  results  are  shown  in  Fig.  2,  in  which  we  have  chosen  the  frequency 
ratio  “/Wp  such  that  there  are  almost  no  surface  plasmons  propagating  on  the 
substrate  surface.  Curves  1  and  2  in  both  cases  (a)  and  (b)  correspond  to  low 
frequencies  (w  -  0.1  and  high  frequencies  (w  -  “p) .  respectively.  Since 
surface  plasmon  waves  are  absent  in  these  cases,  the  enhancement  of  the  decay 
rate  in  every  case  shown  in  the  figures  must  be  due  to  absorption  by  metallic 
granules.  We  attempt  below  to  understand  qualitatively  the  peaks  in  these 
curves . 

In  the  low-frequency  regime,  the  absorption  peak  appears  at  around  the 

conduction  percolation  threshold  q^  of  metallic  granules,  qualitatively  in 

agreement  with  what  has  been  found  in  Ref.  3  in  which  a  simple  model  of  q^  - 

1/3  is  assumed.  This  is  understandable  because  only  at  q  -  q^  there  is 

sufficient  metal  component  to  absorb  the  energy  transferred  from  the  excited 

admolecule.  As  q  increases  further,  metallic  particles  become  closer  to  a 

bulk  metal  and  the  absorption  rate  is  reduced  again.  This  is  because  the 

collisions  of  carriers  responsible  for  the  absorption  are  usually  more 

frequent  in  metallic  particles  than  in  the  bulk.  In  the  high-frequency 

regime,  the  position  of  the  absorption  peak  can  be  very  different  for 

different  q  .  When  q  -  0.3,  curve  2  peaks  at  about  q  -  0.7  in  Fig.  1(a),  and 
c  c 

when  q  -  0.5,  it  peaks  at  q  in  Fig.  1(b).  To  account  for  this  interesting 

c  c 
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phenomenon,  we  believe  that  in  the  optical  wavelength  range  there  exist  two 

distinct  percolation  thresholds,  where  one  results  from  the  usual  percolation 

formed  by  one  connected  path  extending  throughout  the  material  and  the  other 

from  the  "surface  percolation"  formed  by  so  many  connected  paths  that  almost 

19 

the  whole  cross  section  is  conducting.  For  low  frequencies,  the  usual 

percolation  dominates  and  the  optical  behavior  of  the  composite  material  is 

2  20 

similar  to  that  of  the  static  case.  ’  For  high  frequencies,  the  bulk  of  the 
metallic  part  in  the  granular  material  is  responsible  for  the  optical 
absorption  in  the  substrate.  Hence  the  peak  of  curve  2  clearly  indicates  the 
existence  of  a  new  percolation  thrfshold  resulting  from  the  bulk  effects. 

When  the  transition  frequency  is  such  that  there  exist  surface  plasmon 
waves  in  the  substrate,  energy  can  be  transferred  from  the  excited  admolecule 
both  by  the  absorption  processes  described  above  and  the  excitation  of  surface 
plasmon  waves.  The  latter  is  found  much  more  important  than  the  former. 
Numerical  results  are  plotted  in  Fig.  3  for  some  typical  cases.  It  is  clearly 
seen  that  the  decay  rate  is  greatly  enhanced,  especially  when  q  ~  0.75  for  q^ 
-  0.3  and  when  q  -  0.85  for  q^  -  0.5.  A  physical  understanding  of  this 
behavior,  however,  is  not  straightforward,  and  further  investigation  is 
required. 

This  research  was  supported  in  part  by  the  Office  of  Naval  Research  and 
the  National  Science  Foundation  under  Grant  CHE-9016789. 
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Figure  Captions 

1.  Effective  dielectric  function  £2^^^)  volume  fraction  of  metallic 

particles  for  -  0.3.  The  solid  (dashed)  lines  represent  the  real 
(imaginary)  part  of  £2(w).  (a)  w/w^  -  0.1,  (b)  -  0.35  and  (c)  -  1.0 

2.  Luminescence  decay  rate  7  of  the  admolecule  as  a  function  of  q.  Curve  1 

is  for  «/Wp  -  0.1  and  curve  2  is  for  -1.0.  (a)  q-0.3  and  (b)  q  - 

0.5. 

3.  Luminescence  decay  rate  7  of  the  admolecule  vs.  q  in  the  frequency  region 

where  the  energy  transfer  from  the  admolecule  to  surface  plasmons  occurs. 

Curve  1:  q  -  0.3,  w/u  -  0.35;  Curve  2:  q  -  0.5,  w/w  =  0.45. 
c  p  c  p 


